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DOSE RESPONSE FUNCTIONS IN THE ATMOSPHERE 
DUE TO INCIDENT HIGH-ENERGY PROTONS WITH APPLICATION 
TO SOLAR PROTON EVENTS 

By John W. Wilson, Jules J. Lambiotte, Jr., 

Trutz Foelsche, and Tassos A. Filippas 
Langley Research Center 

SUMMARY 

The results of a recent high-energy nucleon-transport calculation are presented 
as dose response functions at five altitudes in the atmosphere. The dose in extremities 
and the dose averaged over the body are directly calculated from these response functions 
by integration over the incident- energy spectra in the range 0.1 to 10 GeV. A dose- 
altitude profile is obtained by using all five sets of response functions. Also presented 
is a geomagnetic- cutoff model for normal magnetic-field strengths. Doses for several 
solar events of solar cycle 19, including the high-energy event of February 23, 1956, are 
presented. The present results are compared with other calculations. 

INTRODUCTION 

Flamm and Lingenfelter (ref. 1) calculated the dose in the atmosphere for idealized 
solar proton spectra. They used rigidity spectra of the form N = Nq exp^-P^Poj, where 
P is particle rigidity, N is the number of protons with rigidity greater than P, and 
Nq and Pq are parameters. ^The parameter Pq assumed values in the range 
75 MV/c to 350 MV/c.) The experimental total nonelastic cross section and multiplicity 
were used. It was assumed in their calculations that the incident energy was equally 
shared among all secondary particles. Their results are of limited usefulness since 
many important solar spectra cannot be well represented by an exponential rigidity spec- 
trum (ref. 2). Furthermore, the energy ranges of the most important solar events are 
not covered by these calculations. 

Leimdorfer et al. (ref. 3) used a rigidity spectrum with Pq = 100 MV/c and a 
Monte Carlo transport code to calculate dose in the atmosphere. They used the model 
of Bertini (ref. 4), which is an intranuclear cascade followed by evaporation, to calculate 
nuclear interactions. Leimdorfer et al. divided the incident- energy spectra from 50 to 
450 MeV into 50 MeV groups and tabulated the dose due to each group. Their tabulations 



are somewhat flexible because a variety of incoming spectra can be treated, but these 
spectra were cut off at an incoming energy of 450 MeV (P < 1000 MV/c). They con- 
sidered the effects of shielding by an additional layer of metal (skin of an airplane) of 
various thicknesses. Their results showed that at atmospheric depths of 58 to 22 g/cm2- 
thin metal shields (1 to 5 g/cm2 of Fe) do not greatly affect the dose. 

Armstrong et al. (ref. 5), using the compilations of Foelsche (as presented in 
ref. 6), have recently completed calculations of dose rates as a function of altitude for 
the upper and lower limits of the February 23, 1956, solar proton event. Their calcula- 
tions include proton energies up to 3 GeV. As will be seen in "Results and Discussion," 
the results presented in reference 5 are higher by a factor of 2 to 3 than the average dose 
calculations presented herein. The calculations of reference 5 used maximum-dose con- 
version factors for normally incident neutrons and skin-dose conversion factors for nor- 
mally incident protons . The two calculations agree within statistical accuracy when the 
same conversion factors are used (unpublished calculations of R. G. Alsmiller, Jr., Oak 
Ridge National Laboratory). 

This paper presents spectrum-independent dose response functions for isotropically 
incident protons in the energy range 0.1 to 10 GeV as a function of altitude. By computing 
dose response rather than dose for a given spectrum, a single transport calculation gives 
all the response functions. These response functions can then be used to calculate dose 
for any number of arbitrary spectra within the energy range 0.1 to 10 GeV. Once the 
dose response functions are available, computations can be done in real time during a 
solar event. Thus, if the event spectrum can be measured, the altitude profile of dose 
rates can be predicted during the event. The shape of the curve of dose response func- 
tion plotted against incident energy also indicates which measurements of the incident 
spectra are needed (that is, which energies are most important) for supersonic transport 
(SST) operations. 

Also presented are the dose rates computed through the application of these 
response functions to several important solar proton events of solar cycle 19. Geomag- 
netic cutoffs have been incorporated by proper modification of the input spectra. 

These response functions have been used previously by Foelsche et al. (ref. 6) to 
calculate dose rates in the upper atmosphere caused by the most important solar proton 
events of the last solar cycle and by galactic cosmic- ray protons. These calculations 
are compared with experimental results for galactic cosmic rays and discussed in 
reference 6. 
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SYMBOLS 



Dt 

E 


photon current-to-dose-rate conversion factor, 

photon/cm2 

speed of light 

dose due to gamma rays (photons), rad 

estimated total dose due to incident high-energy protons, rem or rad 
energy, GeV 


E 


y 


Ev 


Ew 

I(E Amag) 


N 

No 

P 

Po 

QF 


photon energy, MeV 

proton cutoff energy from vertical, GeV 

proton cutoff energy from western horizon, GeV 

photon flux density, photons/cm2-sec 

transmission coefficient for protons of energy E at magnetic latitude 
Xmag to Earth's magnetic field, dimensionless 

idealized exponential rigidity flux density of solar protons, protons/cm2-hr 

total flux density of solar protons, protons/cm2-hr 

proton rigidity, MV/c 

rigidity parameter, MV/c 

quality factor, rem/rad 


Rn(E) dose response due to secondary neutrons produced by protons of energy 

rem or rad 


incident on top of atmosphere. 


incident proton/cm2 
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Rp(E) dose response due to primary and secondary protons from protons of 

rem or rad 


energy E incident on top of atmosphere, 


incident proton/cm2 


Rt(E) 


^mag 

$(E) 


total dose response function due to protons of energy E incident on 

top of atmosphere, 

incident proton/cm2 

photon detector efficiency, counts/photon 
magnetic latitude, deg 

proton fluence (through plane) incident on top of atmosphere, protons /cm2-GeV 


cp^(E^ differential-energy flux density, photons/cm2-sec-MeV 

n solid angle, sr 


METHOD OF CALCULATION 
The Transport Code 

The transport code is a set of computer programs written for the CDC 6600 com- 
puter. These programs record the history of each incident particle and its progeny until 
they are stopped, absorbed, or thermalized. An analysis program then reads the history 
tapes and compiles statistics on the fate of each generation of particles. For example, 
the statistics compiled for this study are differential- energy flux densities for neutrons 
and protons at various altitudes . 

The transport program for energies below 400 MeV was written by Leimdorfer 
et al. and is described in reference 7. The basic structure of the Langley Research 
Center program is the same as that in reference 7. The Langley program is an exten- 
sion to the GeV range and includes the transport of pions. This extension required 
nuclear-interaction data which were obtained as described in the following sections. 

Range 0.75 to 2 Ge V.- The nuclear- interaction data in this range were calculated 
by Bertini (ref. 8). Bertini's calculations of nuclear cross sections employ an intranu- 
clear cascade in an impulse approximation. The impulse approximation assumes that 
the nucleus at all times during the interaction is composed only of nucleons (i.e., no 
exchange currents or resonances). The incident nucleons and resultant cascade nucleons 
are furthermore assumed to have only binary, or two-body, interactions, where the dis- 
tributions of nucleons within the nucleus are taken from electromagnetic form factors. 
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Bertini assumed that pion production occurred through the 



resonance, where the 


decay, in the center- of-momentum system, was 50 percent isotropic, 25 percent forward, 
and 25 percent backward. 


The proton- interaction data were taken from reference 8 . This reference contains 
data for protons of energies 0.75, 1, and 2 GeV, and for neutrons of energy 1 GeV on the 
target elements and Pb207. 

The neutron- interaction data within the same energy range were not available. The 
neutron data at energies of 0.75, 1, and 2 GeV on 0^6 were calculated for the present 
paper by using SU(2) symmetry and the proton data of reference 8 . (Note that has 
zero isotopic spin and is therefore a SU(2) scalar.) The cross sections, multiplicities, 
and energy-angle distributions found for neutrons (by using SU(2)) agreed to within 
10 percent with the neutron data of reference 8 at 1 GeV. Thus, the interaction data for 
neutrons to 2 GeV were completed, and the Bertini calculations were shown to be consis- 
tent with this basic symmetry principle of strong interactions. 


Range 2.0 to 10 GeV.- The range of energy for which nuclear interaction data are 
known (to 2 GeV) is not sufficient for solar events such as that which occurred on 
February 23, 1956, or especially for galactic cosmic rays. Since the experimental 
specific-yield functions for low-energy neutrons indicate that the energy dependence of 
primary protons in the GeV range only slightly affects the secondary nuclear yields in 
the upper atmosphere (ref. 9), the multiplicities and energy-angle distributions (normal- 
ized to the incident energy) were assumed not to change from 2 to 10 GeV. The pion 
component will, of course, not be correct at the higher energies since the onset of 
diffraction-production is ignored. Because energy is absorbed when pions are produced, 
the secondary nucleons will have a little less energy than the calculations predict. The 
pions are assumed not to have nuclear interactions since at high altitudes, their mean 
free path to decay is much smaller than their mean free path to nuclear interaction. 

The transport calculations were carried out for 12 500 protons uniformly distributed 
over 25 incident-energy groups to cover the range 0.1 to 10 GeV. The final particle flux 
densities were compiled from approximately four million particle events, including those 
as many as 12 generations removed from the incident primary protons. 


The Atmospheric Model 

The thickness of the traversed layer in the atmosphere is expressed in units of 
g/cm2. In this way the need for expressing explicitly the change of atmospheric density 
with depth is eliminated. Therefore, the model used for the atmosphere was considered 
to be an infinite layer of air of uniform density. The radiation was assumed to be iso- 
tropic (see ref. 10 for discussion) from the upper hemisphere. 
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The perturbation of the radiation field due to the presence of tissue and the airplane 
was neglected in the transport calculations. The length of the airplane, approximately 
100 m, is small compared with the mean free path, which is greater than 1 km at SST 
cruising altitudes. Hence, because of the isotropy of the incident radiation and the large 
mean free path of protons and neutrons at even moderate energies, the penumbra formed 
will completely overcast the shadow of the entire airplane. 

Dose Response Functions 

Current-to-dose-rate conversion factors were used to transcribe flux densities 
obtained from the transport calculations to dose response. The conversion factors used 
were for a tissue composition of H^, n 14^ and Below 60 MeV, the factors for 

protons were calculated by assuming no nuclear interactions. The data for protons from 
60 to 400 MeV were taken from Turner et al. (ref. 11). A simple energy- loss calculation 
indicates that the absorbed-dose conversion factors above 400 MeV can be assumed con- 
stant and equal to the absorbed-dose conversion factor of reference 11 at 400 MeV. The 
rem (dose equivalent) conversion above 400 MeV was found by using the average quality 
factor (ratio of rem to rad, see ref. 12) of 1.4 of Turner et al., which was computed at 
400 MeV with nuclear interactions. The neutron current-to-dose-rate conversion factors 
were taken from Irving et al. (ref. 13) and Kinney and Zerby (ref. 14). All these conver- 
sion factors include the nuclear-star dose equivalent (nuclear interactions in tissue with 
multiple low-energy prongs) from protons and neutrons. The average- and skin-dose 
conversion factors for a 30-cm-thick slab of tissue for isotropic neutrons and protons 
were used in the present calculations. The extremity dose response functions were 
obtained by using the omnidirectional flux density and skin-dose conversion factors. 

(See ref. 6 for discussion.) 

The contribution to the total dose due to particles other than nucleons (e.g., pions, 
electrons, and gammas) has been neglected. The absorbed dose (rad) will be too low by 
possibly as much as 30 percent but the dose equivalent (rem) will not be so greatly affected 
because of the low biological importance of these neglected radiations (that is, the quality 
factor QF ~ 1). For example, the dose from y-rays produced by galactic cosmic rays, 
estimated from the measurements of Haymes (ref. 15) at the y-transition maximum, is 
computed from 

where is the photon current-to-dose-rate conversion factor and (p^ is the photon 

differential-energy flux density. In magnitude 
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C < 10-6 mrad 
^ photon/cm2 


Define 


p3.5 MeV 
•^0.5 MeV 



where fy is known from reference 15 except for detector efficiency rjy Thus from 
reference 15, the photon flux density in the energy range 0.5 to 3.5 MeV is 


f = 1.8 X 1Q3 photons 
^ Vy cm2-hr 

and the dose is overestimated by 


D 


y 


~ 10-6 


/l.8 X 103^ 

mrad 

1 ^y y 

hr 


fl.8 X 1Q-3\ 

V / 


mrad 

hr 


If 


'y 


is on the order of 10 percent, then 


Dy § 0.018 

^ hr 


This is on the order of 3 percent of the total absorbed dose, which was 0.65 mrad/hr. 
The dose due to y-rays D.^ is only about 1 percent of the total dose equivalent 
(1.3 mrem/hr) at the same altitude. (See ref. 6.) Note that with = 10 percent, the 
total dose from y-rays is on the order of 10 percent of the neutron absorbed dose rate 
(0.2 mrad/hr) of reference 6. 


DOSE CALCULATION 


Geomagnetic Cutoff 

The normal geomagnetic- cutoff model can be represented by 


l(E,Xmag) = 1 - exp(E^ 

A(Amag) = 0-693 

V (Ey - E^) 


- EjA^Xmag) 


where l(E,Xjnag) fraction of intensity transmitted at energy E, and Ey and Ey; 

are the vertical and western cutoff energies, respectively, given in figure 1 (data taken 
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from ref. 16). This model does not account for the change in cutoff due to depression of 
the earth's magnetic field by solar cosmic rays, a situation which oftens occurs during 
periods of high solar activity. This change in cutoff has not been included since it is as 
yet a not too well understood phenomenon. 

Dose Integral 

With the dose response functions, the dose averaged over the body or the dose in 
extremities is easily calculated when the incident spectra are known. The dose at high 
latitudes A^ag > is given by 

plO GeV 

Dt = \ Rt(E) $(E)dE 

-"O.! GeV 

where Rt(E) is the corresponding dose response function and 4> is the incident proton 
fluence through a plane* in protons/cm2-GeV. To include latitude effects, then 

plO GeV 

^t^^mag) “ ^ I^t(E)I^E,X);nag)^(E)dE 

where l(E,Xniag) the transmission at geomagnetic latitude Xmag and was given in 
the previous section. 


RESULTS AND DISCUSSION 

The dose response functions for dose averaged over the body and for extremities 
are shown for altitudes from 300 to 20 g/cm2 in figures 2 to 7. These curves were 
obtained by smoothing the Monte Carlo data given in tables I to VI. The greatest statisti- 
cal fluctuations were at the lower incident energies of the proton dose response functions 
(figs. 2 and 5), where the curves are falling by orders of magnitude. The neutron dose 
response had, in general, considerably less statistical scatter than that of protons. This 
is due to the much larger number of neutrons, which provides a larger statistical sample. 
The slight depression in dose at 0.4 GeV as seen in the tables is due to the fact that two 
different transport codes joined the data at 0.4 GeV. These codes usually give slightly 
different results at the energy where they are combined. 

At altitudes for subsonic flight (about 300 g/cm2) no significant dose is generally 
expected. At the present SST cruising altitudes of 19.8 km (about 58 g/cm2), moderate 

* 1 

Fluence through a plane = - fluence through a sphere: 

$ protons/cm2-GeV = ^ 2tt sr ^ protons/cm2-GeV-sr. 

<5 dii 
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radiation doses are expected from solar events with energies above 0.25 GeV. At 
cruising altitudes of 20 g/cm2, virtually all solar events (except the low-energy events 
E ~ 0.1 GeV) have potential for producing a significant dose. This is particularly true 
for the extremity dose at 20 g/cm2 as seen in figures 5 and 7. The peak in the dose 
response function (Bragg peak) is due to uncollided primary protons. This peak does 
not occur at the lower altitudes shown in the figures since nearly all primary protons 
have suffered at least one nuclear interaction. 

The surface dose as calculated in reference 1 is presented in figure 8 in compari- 
son with the extremity dose of the present work for rigidities of 150, 250, and 350 MV/c. 
The greatest differences at lower altitudes are presumably due to the simplified nuclear 
interaction model used in reference 1; in this model the incident energy at interaction 
was assumed to be equally shared among all secondary particles. 

The presently calculated average doses are compared in figure 9 with the results 
of reference 3 for a rigidity spectrum with Pq = 100 MV/c. The absorbed dose (rad) of 
reference 3 is in good agreement with the present calculations; however, at lower altitudes 
the dose equivalent (rem) is different. This indicates some discrepancies in the calcula- 
tion of dose for neutrons in the two calculations. 

Spectra for two solar events of cycle 19 are shown in figure 10 as compiled by 
Foelsche from measured data (ref. 6). The prompt spectra of the February 23, 1956, 
event were observed on different locations during the maximum phase of the event in 
the first hours after particle onset. The spectra of the November 12, 1960, event were 
observed 5 hours (1840 UT), 10 hours (2330 UT), and 27 hours (1603 UT on Nov. 13) after 
the particle onset. Extensive results (for whole-body average dose rates) of the present 
calculations applied to these spectra can be found in reference 6. 

The extremity dose rates for the upper and lower limits of the February 23, 1956, 
event are shown in figure 11 along with similar calculations from reference 5. In refer- 
ence 5, the spectra were cut off above 3 GeV and the conversion factors used were skin 
dose for protons and maximum dose for neutrons. Although the transport calculations 
were for protons isotropically incident at the top of the atmosphere, the conversion fac- 
tors used in reference 5 were for neutrons and protons normally incident on the tissue 
sample. The agreement of the two results is reasonably good, especially for the lower 
limit to the February 1956 event. The average dose rate for the same event as presented 
in reference 6 is lower by a factor of 2 to 3 than the extremity dose rates of the present 
calculations and the dose rates of reference 5. When the same conversion factors are 
used and the primary spectra are cut off at 3 GeV, the calculations of reference 5 and the 
present calculations are in complete agreement. The average and extremity dose rates 
during the November 12, 1960, event are shown in figures 12 and 13. This event is char- 
acteristic of medium-energy, high- intensity solar events. 
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CONCLUDING REMARKS 


The importance of the dose response functions with respect to the supersonic trans- 
port (SST) is twofold. First, they have been used by Foelsche et al. to assess theoreti- 
cally the radiation hazard to the SST. There is, of course, the need for experimental 
verification of these results. Second, with the response functions presented herein and 
with sufficient real-time measurements by both satellite and by ground-based neutron 
monitors, it is possible that current estimates of extremity and depth dose rates can be 
given so that SST flights may be rescheduled, or evasive measures can be taken to mini- 
mize exposure. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 2, 1970. 
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TABLE I.- PROTON AVERAGE DOSE RESPONSE FUNCTIONS 
DUE TO INCIDENT PROTONS 


Energy 

range. 

20 g/cm2 

58 g/cm2 1 

100 g/cm2 I 

200 g/cm2 j 

300 g/cm2 




rad 


GeV 



incident proton/cm's 

.10-. 15 

0. 

0. 

0. 

0. 

0. 

.15-. 20 

2.75E-09 

0. 

0. 

0. 

0. 

.20-. 25 

1 . 78E-08 

5. 19F-11 

0. 

0. 

0. 

.25-. 30 

3.38E-08 

1.04E-10 

0. 

0. 

0. 

.30-. 35 

5.21E-08 

2.01E-09 

0. 

0. 

0. 

.35-. 40 

6.05E-08 

1. 16E-08 

1.04E-10 

0. 

0 . 

.40-. 50 

4.02E-08 

1.37E-08 

1.78E-09 

0. 

0. 

,50-. 60 

7.88E-08 

3.68E-08 

1. 16E-08 

0. 

0. 

.60-. 70 

8.80E-08 

4.78E-08 

2.32E-08 

4.936-10 

0. 

,70-. 80 

8.76E-08 

4. 86E-08 

2.91E-08 

3.046-09 

1 .956-10 

.80-. 90 

9.35E-08 

5. 676-08 

2.99E-08 

6.306-09 

4.816-10 

.90-1.0 

9.24E-08 

5.686-08 

3.716-08 

1 .026-08 

1 .026-09 

1.0-1. 5 

9.94E-08 

6. 72E-08 

4.32E-08 

1.706-08 

5.216-09 

1. 5-2.0 

1 .09E-07 

8.33E-08 

6.03E-08 

2.706-08 

1 .016-08 

2.0-2. 5 

1 .12E-07 

8.32E-08 

6.326-08 

2 . 766-08 

1.246-08 

2. 5-3.0 

1. 19E-07 

9. 35E-08 

7.236-08 

3.766-08 

1 .676-08 

3.0-3. 5 

1 .27E-07 

1.13E-07 

8.33E-Q8 

4.046-08 

2.236-08 

3. 5-4.0 

1 .40E-07 

1. 076-07 

9.22E-08 

5.496-08 

2.816-08 

4. 0-4. 5 

1 .45E-07 

1.326-07 

1.056-07 

6.226-08 

2.836-08 

4. 5-5.0 

1 .62E-07 

1.476-07 

1.156-07 

6.096-08 

3.356-08 

5. 0-6.0 

1 .68E-07 

1.516-07 

1.236-07 

7.586-08 

3.716-08 

6. 0-7.0 

1 .88E-07 

1. 72E-07 

1.396-07 

9.196-08 

4.956-08 

7. 0-8.0 

1 .94E-07 

1.866-07 

1.716-07 

1.056-07 

5.846-08 

8. 0-9.0 

2.10E-07 

2. 02E-07 

1.626-07 

1.256-07 

6.886-08 

9.0-10.0 

2 .41E-07 

2.31E-07 

2. 056-07 

1.176-07 

7.266-08 



rem 



incident proton/cm4 

.10-. 15 

0. 

0. 

0. 

0. 

0. 

.15-. 20 

3 .99E-09 

0. 

0. 

0. 

0. 

.20-. 25 

2 .54E-08 

7.52E-1 1 

0. 

0. 

0. 

.25-. 30 

4 .69E-08 

1.506-10 

0. 

0. 

0. 

.30-. 35 

6.90E-08 

2.896-09 

0. 

0. 

0. 

.35-. 40 

8. lOE-08 

1.646-08 

1.506-10 

0. 

0. 

.40-. 50 

5.59E-08 

1.846-08 

2. 536-09 

0. 

0. 

.50-. 60 

1 . llE-07 

5.096-08 

1.596-08 

0. 

0. 

,60-. 70 

1.24E-07 

6.686-08 

3.226-08 

7.066-10 

0. 

.70-. 80 

1 .24E-07 

6. 86E-08 

4.076-08 

4.206-09 

2 .786-10 

.80-. 90 

1 .33E-07 

8. OlE-08 

4. 196-08 

8.666-09 

6.836-10 

.90-1.0 

1 .31E-07 

8.02E-08 

5.216-08 

1 .436-08 

1 .406-09 

1.0-1. 5 

1 .416-07 

9.46E-08 

6. 076-08 

2.366-08 

7.316-09 

1 .5-2.0 

1 .55E-07 

1 .186-07 

8. 506-08 

3.816-08 

1 .386-08 

2.0-2. 5 

1 .58F-07 

1.1 76-07 

8.916-08 

3.876-08 

1.746-08 

2. 5-3.0 

1 .68E-07 

1. 326-07 

1.026-07 

5.286-08 

2.336-08 

3.0-3. 5 

1.80E-07 

1.59E-07 

1. 176-07 

5.676-08 

3.146-08 

3. 5-4.0 

1 .98E-07 

1.51E-07 

1.306-07 

7.726-08 

3.926-08 

4.0-4. 5 

2.06E-07 

1.86E-07 

1.486-07 

8.726-08 

3.936-08 

4. 5-5.0 

2.28E-07 

2.07E-07 

1.626-07 

8.586-08 

4.696-08 

5. 0^6.0 

2.38E-07 

2.1 26-07 

1.736-07 

1.076-07 

5.176-08 

6. 0-7.0 

2.66E-07 

2.42E-07 

1.966-07 

1.296-07 

6.916-08 

7. 0-8.0 

2.74E-07 

2. 626-07 

2.406-07 

1.476-07 

8.186-08 

8. 0-9.0 

2 .96E-07 

2.84E-07 

2. 296-07 

1 .756-07 

9.666-08 

9.0-10.0 

3.40E-07 

3.256-07 

2.886-07 

1 .656-07 

1 ,026-07 



TABLE n.- NEUTRON AVERAGE DOSE RESPONSE FUNCTIONS 
DUE TO INCIDENT PROTONS 


Energy 

20 g/cm2 

58 g/cm2 

100 g/cm2 

200 g/cm2 

300 g/cm2 

range, 




rad 


GeV 

Dose response function, 

— 




incident proton/cm^ 

. 10- . 1 5 

6.21E-10 

3. 835-10 

1. 705-10 

0. 

0. 

.15-. ?0 

1 .OOE-09 

4.405-10 

1.995-10 

5.005-11 

0. 

. 20-.25 

1 .64E-09 

1.025-09 

6.965- 10 

2.445-10 

8.625-11 

.25-. 30 

2.315-09 

1.405-09 

7.835-10 

2.335-10 

1 .055-10 

.30-. 35 

2.775-09 

2.455-09 

1. 51 5-09 

3.895-10 

1 .635-10 

.35-.^0 

3 . 745-09 

3.165-09 

2.075-09 

6.135-10 

1.635-10 

.A0-.50 

3.245-09 

3. 1 15-09 

2.245-09 

7.705-10 

2.635-10 

.50-. 60 

8. 365-09 

8.615-99 

6.465-09 

2.715-09 

1 .145-09 

.60-. 70 

9.585-09 

1.105-08 

8.615-09 

4.005-09 

1 .655-09 

. 70-.80 

1 . 225-08 

1.455-08 

1.215-08 

5.755-09 

2.465-09 

. 80- . 90 

1 .275-08 

1.515-08 

1.465-08 

7.805-09 

3.815-09 

.90-1.0 

1.465-08 

1. 595-08 

1.575-08 

9.935-09 

4.955-09 

1 .0-1.5 

1. 815-08 

2. 165-08 

1.995-08 

1.275-08 

7.155-09 

1 .5-2.0 

2 .255-08 

2.565-08 

2.515-08 

1 .795-08 

1.175-08 

2. 0-2. 5 

2 . 735-08 

3. 255-08 

3.205-08 

2.415-08 

1.345-08 

2. 5-3.0 

3 .035-08 

3. 735-08 

3.675-08 

2.835-08 

1.785-08 

3 .0-3.5 

3.82F-08 

4.495-08 

4.395-08 

3.285-08 

2.255-08 

3. 5-4.0 

3.815-08 

4.575-08 

4.635-08 

3.765-08 

2 .625-08 

4. 0-4. 5 

4 . 1 85-08 

5. 425-08 

5.455-08 

4.325-08 

2.925-08 

4. 5-5.0 

4.735-08 

5.985-08 

5.895-08 

4.685-08 

2.945-08 

5. 0-6.0 

5.425-08 

6. 515-08 

6.665-08 

5.435-08 

3.655-08 

6. 0-7.0 

6.225-08 

7.825-08 

7.915-08 

6.425-08 

4.565-08 

7. 0-8.0 

7 .105-08 

8. 945-08 

9.085-08 

7.285-08 

4.995-08 

a. 0-9.0 

8 .165-08 

9. 695-08 

1.025-07 

8.445-08 

6.015-08 

9.0-10.0 

8. 775-08 

1. 105-07 

1 .135-07 

9.265-08 

6.175-08 


Dose response function. 

rem 





incident proton/cm2 

. I0-. 15 

4 . 005-09 

2.455-09 

1.105-09 

0. 

0. 

.15-. 20 

6.1 15-09 

2.595-00 

1.075-09 

2.795-10 

0. 

. 20-.25 

9.325-09 

5. 745-09 

3.795-09 

1.295-09 

4.535-10 

.25-. 30 

1 .2 75-08 

7. 755-09 

4.385-09 

1.295-09 

6.125-10 

.30-. 35 

1 .475-08 

1.305-08 

8.245-09 

2.315-09 

9.725-10 

.35-. 40 

1 .955-08 

1.675-08 

1.065-08 

3.185-09 

6.855-1 0 

.40-. 50 

1 . 505-08 

1.565-08 

1.145-08 

3.725-09 

1.215-09 

. 50- . 60 

4.00E-Q8 

4.045-08 

3.045-08 

1.335-08 

5.685-09 

.60-. 70 

4.715-08 

5.195-08 

4.005-08 

1.995-08 

8.395-09 

.70- . 30 

6.025-08 

7.015-08 

5.985-08 

2.645-08 

1.165-08 

.80-. 90 

5 .355-08 

7. 265-08 

6.795-08 

3.685-08 

1.895-08 

.90-1.0 

7.415-08 

7. 735-08 

7. 395-08 

4.665-08 

2.265-08 

1.0-1. 5 

9.375-08 

1.065-07 

9.665-08 

6.025-08 

3.465-08 

1. 5-2.0 

1 .195-07 

1. 315-07 

1.265-07 

8.515-08 

5.675-08 

2. 0-2. 5 

1 .41E-C7 

1.655-07 

1.595-07 

1 .165-07 

6.395-08 

2. 5-3.0 

1 .585-07 

1.905-07 

1.855-07 

1.365-07 

8.645-08 

3. 0-3. 5 

2.035-07 

2. 285-07 

2.195-07 

1.615-07 

1.105-07 

3. 5-4.0 

1 . 955-07 

2.325-97 

2.305-07 

1.835-07 

1 .235-07 

4.0-4. 5 

2.195-07 

2.735-97 

2.695-07 

2.095-07 

1.425-07 

4. 5-5.0 

2.415-07 

2.995-97 

2.925-07 

2.335-07 

1 .455-07 

5. 0-6.0 

2.785-07 

3.315-07 

3.355-07 

2.655-07 

1 . 765-07 

6. 0-7.0 

3.175-07 

3.975-07 

3. 94 5-07 

3.115-07 

2.175-07 

7. 0-8.0 

3 .6 15-07 

4.525-07 

4.515-07 

3.575-07 

2.445-07 

8 .0-9.0 

4.155-07 

4. 865-07 

5.095-07 

4.045-07 

2.875-07 

9. 0-10. 0 

4.455-07 

5.525-07 

5.565-07 

4.515-07 

2.965-07 
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II 



TABLE m.- TOTAL AVERAGE DOSE RESPONSE FUNCTIONS 
DUE TO INCIDENT PROTONS 


Energy 

20 g/cm2 

58 g/cm2 

100 g/cm2 

200 g/cm2 

300 g/cm2 

range, 

GeV 

Dose response function. 

rad 





incident proton/cm2 

.10-. 15 

6.21F-10 

3.83E-10 

1.70E-10 

0. 

0. 

.15-. 20 

3 .75E-09 

4.40E-I0 

1.99E-10 

5.00E-11 

0. 

.20-. 25 

1 .95E-08 

1.07E-09 

6.96E-10 

2.44E-10 

8 .62E-1 1 

.25-. 30 

3.61E-08 

1. 50E-09 

7.83E-10 

2.33E-10 

1.05E-10 

.30-. 35 

5 .49E-08 

4.46E-09 

1.51E-09 

3.89E-10 

1 .63E-10 

.35-. 40 

6.43E-08 

1.48E-08 

2.18E-09 

6.13E-10 

1 .63E-10 

.40-. 50 

<i .34E-08 

1. 69E-08 

4.02E-09 

7.70E-10 

2.63E-10 

.50-. 60 

8.72E-08 

4.54E-08 

1.81E-08 

2.71E-09 

1 .14E-09 

.50-. 70 

9.75E-08 

5.88F-08 

3. 18E-08 

4.49E-09 

1 .65E-09 

.70-. 80 

9 .98E-08 

6.32E-08 

4.12E-08 

8.79E-09 

2.65E-09 

.80-. 90 

1.06E-07 

7. 19E-08 

4.45E-08 

1 .41E-08 

4 .29E-09 

.90-1.0 

1.07E-07 

7. 27E-08 

5.28E-08 

2.02E-08 

5.97E-09 

1 .0-1.5 

1.18E-07 

8. 88E-08 

6.31E-08 

2.96E-08 

1 .24E-08 

1. 5-2.0 

1.32E-07 

1.09E-07 

8.54E-08 

4.50E-08 

2.17E-08 

2. 0-2. 5 

1 .39E-07 

1. 16E-07 

9.52E-08 

5.16E-08 

2.58E-08 

2. 5-3.0 

1 .49E-07 

1. 31E-07 

1.09E-07 

6.59E-08 

3 .45E-08 

3. 0-3. 5 

1 .66E-07 

1. 58E-07 

1.27E-07 

7.32E-08 

4.48E-08 

3. 5-4.0 

1 . 78E-07 

1. 53E-07 

1.38E-07 

9.25E-08 

5.43E-08 

4.0-4. 5 

1 .87E-07 

1 .86E-07 

1.60E-07 

1 .05E-07 

5.75E-08 

4. 5-5.0 

2 .09E-07 

2.07E-07 

1.74E-07 

1 .08E-07 

6.29E-08 

5. 0-6.0 

2 .22E-07 

2. 16E-07 

1.90E-07 

1 .30E-07 

7.36E-08 

6.0-7. 0 

2.50E-07 

2. 50E-07 

2.18E-07 

1.56E-07 

9.51E-08 

7. 0-8.0 

2.65E-07 

2. 75E-07 

2.62E-07 

1 .77E-07 

1.08E-07 

8. 0-9.0 

2.91E-07 

2. 99E-07 

2.64E-07 

2.09E-07 

1 .29E-07 

9.0-10.0 

3.29E-07 

3.41E-07 

3.17E-07 

2.10E-07 

1 .34E-07 



rem 





' incident proton/cm2 

. 10-. 15 

4. OOE-09 

2.45E-09 

l.lOE-09 

0. 

0. 

. 15-.20 

1 .OlE-08 

2.59E-09 

1.07E-09 

2.79E-10 

0. 

.20-. 25 

3 .47E-08 

5.81E-09 

3.79E-09 

1 .29E-09 

4.53E-10 

.25-. 30 

5 .95E-08 

7.90E-09 

4.38E-09 

1 .29E-09 

6.12E-10 

.30-. 35 

8.37E-08 

1.59E-08 

8.24E-09 

2.31 E-09 

9.72E-I0 

.35-. 40 

1 .OOE-07 

3.30E-08 

1.08E-08 

3.18E-09 

6.85E-10 

.40-. 50 

7.09E-08 

3.40E-08 

1.40E-08 

3.72E-09 

1 .21E-09 

.50-. 60 

1.51E-07 

9.14E-08 

4.62E-08 

1.33E-08 

5.68E-09 

.50-. 70 

1.72E-07 

1. 19E-07 

7.22E-08 

2.06E-08 

8.39E-09 

. 70- . 80 

1 .84E-07 

1.39F-07 

l.OOE-07 

3.06E-08 

1.19E-08 

.80-. 90 

1 .96E-07 

1.53E-07 

l.lOE-07 

4.54E-08 

1 .95E-08 

.90-1.0 

2.05E-07 

1. 58E-07 

1.26E-07 

6.08E-08 

2.40E-08 

1.0-1. 5 

2.35E-07 

2.01E-07 

1.57E-07 

8.38E-08 

4.19E-08 

1. 5-2.0 

2.74E-07 

2.49E-07 

2.11E-07 

1 .23E-07 

7.04E-08 

2.0-2. 5 

2 .99E-07 

2. 83E-07 

2.48E-07 

1 .55E-07 

8.12E-08 

2. 5-3.0 

3.26E-07 

3. 22F-07 

2.87E-07 

1 .89E-07 

l.lOE-07 

3.0-3. 5 

3.83E-07 

3.87E-07 

3.36E-07 

2.18E-07 

1.41E-07 

3. 5-4.0 

3 .93E-07 

3.83E-07 

3.61E-07 

2.50E-07 

1 .62E-07 

4.0-4. 5 

4.25E-07 

4.59E-07 

4.17E-07 

2.95E-07 

1 .82E-07 

4. 5-5.0 

4.69E-07 

5.05E-07 

4.54E-07 

3.18E-07 

1.92E-07 

5. 0-6.0 

5.16E-07 

5.43E-07 

5.08E-07 

3.71E-07 

2.28E-07 

5. 0-7.0 

5.83E-07 

6.39E-07 

5.90E-07 

4.40E-07 

2.86E-07 

7. 0-8.0 

6.35E-07 

7.14E-07 

6.91E-07 

5.04E-07 

3.25E-07 

8. 0-9.0 

7.11E-07 

7.71E-07 

7.38E-07 

5.80E-07 

3.83E-07 

9.0-10.0 

7.86E-07 

8. 77E-07 

8.44E-07 

6.16E-07 

3.98E-07 



TABLE IV.- PROTON EXTREMITY DOSE RESPONSE FUNCTIONS 
DUE TO INCIDENT PROTONS 


20 g/cm2 



58 g/cm2 


100 g/cm2 


Dose response function, 7 


200 g/cm2 

rad 


300 g/cm2 


incident proton/cm2 


. 10-. 15 

0. 

0. 

0. 

0. 

0. 

. 15-.20 

9.406-08 

0. 

0. 

0. 

0. 

. 20-.25 

1 .546-07 

1.776-09 

0. 

0. 

0. 

.25-. 30 

1.466-07 

3.556-09 

0. 

0. 

0. 

.30-. 35 

1.386-07 

4.646-08 

0. 

0. 

0. 

.35-. 40 

1.206-07 

7.876-08 

3.556-09 

0. 

0. 

.40-. 50 

5.846-08 

4.296-08 

1.166-08 

0. 

0. 

.50-. 60 

1.466-07 

9.136-08 

5. 166-08 

0. 

0. 

. 60-. 70 

1 . 566-07 

1.066-07 

5.706-08 

7.926-09 

0. 

.70-.ro 

1 .656-07 

1.086-07 

6.666-08 

1.146-08 

2. 

.80-. 90 

1 .666-07 

1.1 16-07 

6.976-08 

2.386-08 

3. 

.90-1 .0 

1.806-07 

1.136-07 

9.266-08 

3.926-08 

7. 

1.0-1. 5 

1 .826-07 

1.396-07 

9.976-08 

4.486-08 

1. 

1. 5-2.0 

2.266-07 

1.876-07 

1.206-07 

6.296-08 

2. 

2. 0-2. 5 

2 . 786-07 

1 .926-07 

1,476-07 

7.526-08 

4. 

2. 5- 3.0 

2 . 786-07 

2.346-07 

1.706-07 

1.016-07 

4. 

3. 0-3. 5 

3.076-07 

2.706-07 

2.366-07 

1 .1 76-07 

7. 

3. 5-4.0 

3.226-07 

2.896-07 

2.586-07 

1.686-07 

1 . 

4. 0-4. 5 

3.906-07 

3.366-07 

2. 696-07 

1.796-07 

6. 

4. 5-5.0 

4.056-07 

3.876-07 

3.126-07 

1.756-07 

1 . 

5. 0-6.0 

4.076-07 

3.956-07 

3.226-07 

2.376-07 

1 . 

6. 0-7.0 

5 .176-07 

4.656-07 

3.896-07 

2.596-07 

1 . 

7. 0-8.0 

5 .196-07 

5.216-07 

4.956-07 

2.966-07 

1 . 

8. 0-9.0 

5.506-07 

6.066-07 

4.626-07 

3.276-07 

1. 

0 

0 

1 

0 

0 

6.506-07 

6.446-07 

5.886-07 

3.286-07 

2. 



rem 









incident proton/cm^ 

.10-. 15 

0. 

0. 

0 . 

0 . 

0. 

.15-. 20 

2.296-07 

0. 

0 . 

0 . 

0. 

.20-. 25 

3.216-07 

4.326-09 

0 . 

0 . 

0. 

.25-. 30 

2.586-07 

8.646-09 

0 , 

0 . 

0. 

.30-. 35 

2 .406-07 

1. 106-07 

0 . 

0 . 

0. 

.35-. 40 

2.036-07 

1.586-07 

8.646-09 

0 . 

0. 

.40-. 50 

1 .166-07 

8.196-08 

2.266-08 

0 . 

0. 

.50-. 60 

2.606-07 

1.686-07 

1.026-07 

0 . 

0. 

.60-. 70 

2.746-07 

1.976-07 

1.006-07 

1.826-08 

0. 

.70-. 80 

2 .956-07 

2.046-07 

1.206-07 

2.046-08 

4. 

.80-. 90 

2.866-07 

1.936-07 

1.306-07 

4.696-08 

5. 

.90-1.0 

3.246-07 

2.036-07 

1.746-07 

7.666-08 

1. 

1 .0-1.5 

3.226-07 

2.536-07 

1.896-07 

8.506-08 

2. 

1. 5-2.0 

4.136-07 

3.466-07 

2.146-07 

1.156-07 

4 . 

2. 0-2. 5 

5.346-07 

3.566-07 

2.726-07 

1.456-07 

9. 

2. 5-3.0 

5.29E-07 

4.466-07 

3.196-07 

1.936-07 

8. 

3. 0-3. 5 

5.836-07 

5. 106-07 

4.636-07 

2.336-07 

1. 

3. 5-4.0 

5.036-07 

5.616-07 

5.026-07 

3.336-07 

2. 

4.0-4. 5 

7.666-07 

6.406-07 

5.136-07 

3.506-07 

1. 

4.5-5. 0 

7.696-07 

7.506-07 

6.046-07 

3.436-07 

2. 

5. 0-6.0 

7.686-07 

7.576-07 

6.116-07 

4.726-07 

2. 

6. 0-7.0 

1 .016-06 

8.996-07 

7.576-07 

4.986-07 

2. 

7 . 0- 8 . 0 

1.016-06 

1.026-06 

9.726-07 

5.766-07 

3. 

8. 0-9.0 

1 .066-06 

1.206-06 

9.076-07 

6.256-07 

3,' 

9.0-10.0 

1.256-06 

1.256-06 

1.146-06 

6. 336-07 

4.^ 


.356-08 


.426-08 

.136-08 

.016-07 

.696-08 


.926-07 
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TABLE V.- NEUTRON EXTREMITY DOSE RESPONSE FUNCTIONS 
DUE TO INCIDENT PROTONS 


Energy 

20 g/cm2 

58 g/cm2 

100 g/cm2 

200 g/cm2 

300 g/cm2 

range, 






GeV 

Dose response function, j 

incident proton/cm' 

? 

.10-. 15 

1 .29E-09 

7.71E-10 

2.81E-10 

0. 

0. 

.15-. 20 

2.72E-09 

1.03E-09 

4.95E-10 

4.86E-11 

0. 

.20-. 25 

3 .91E-09 

2.36E-09 

1. 16E-09 

4.50E-10 

1 .90E-10 

.25-. 30 

5.94E-09 

3. 14E-09 

1. 54E-09 

6.12E-10 

2.35E-10 

.30-. 35 

6.1 3E-09 

4. 57E-09 

2.36E-09 

9.89E-10 

4.23E-10 

.35-. 40 

8.52E-09 

6. 17E-09 

3.36E-09 

1.21E-09 

1.97E- 10 

.40-. 50 

5 .66E-09 

6.14E-09 

4.25E-09 

1.30E-09 

3.62E-10 

.50-. 60 

1 .70E-08 

1.61E-08 

1.13E-08 

3.84E-09 

1 .97E-09 

.60-. 70 

1 .93E-08 

2.06E-08 

1.49E-08 

5.42E-09 

2.60E-09 

.70-. 80 

2 .52F-08 

2. 72F-08 

2.31E-08 

9.91E-09 

3.69F-09 

.80- .90 

2 .83E-08 

2.98E-08 

2.58F-08 

1.27E-08 

6.11E-09 

.90-1.0 

3.29E-08 

3. 32E-08 

2.91E-08 

1 .59E-08 

7.63E-09 

1.0-1. 5 

4.09E-08 

4. 59E-08 

3.98E-08 

2.09E-08 

1.13E-08 

1 .5-2.0 

5.57E-08 

5.90E-08 

5.42E-08 

3.05E-08 

2.09E-08 

2. 0-2. 5 

6 .30E-08 

7. 1 8E-08 

6. 51E-08 

4.40E-08 

2.25E-08 

2. 5-3.0 

7.06E-08 

8. 54E-08 

7.69E-08 

5.14E-08 

3.19E-08 

3.0-3. 5 

8.64E-08 

1 . 03E-07 

9.45E-08 

6.23E-08 

3.98E-08 

3. 5-4.0 

8.47E-08 

1. OlE-07 

l.OOE-07 

7. lOE-08 

4.16E-08 

4. 0-4. 5 

9.73E-08 

1.21E-07 

1.14E-07 

8.15E-08 

5.29F-08 

4. 5-5.0 

1 .C4E-07 

1.32E-07 

1.30E-07 

8.97E-08 

5.29E-08 

5. 0-6.0 

1.20E-07 

1.46F-07 

1.41E-07 

1.03E-07 

6.62E-08 

6. 0-7.0 

1 .32E-07 

1. 72F-07 

1.67E-07 

1 .16E-07 

7.69E-08 

7. 0-8.0 

1.53E-07 

1. 99E-07 

1.85E-07 

1 .35E-07 

9.40E-08 

8. 0-9.0 

1.81E-07 

2.15E-07 

2.14E-07 

1.57E-07 

l.OlE-07 

9.0-10.0 

1.88E-07 

2.35E-07 

2.33E-07 

1.77E-07 

1.12F-07 


Dose response function. 

rem 





incident proton/cm2 

.10-. 15 

6.'*5E-09 

3.87E-09 

i.'iiE-og 

0. 

0 . 

.15-. 20 

1 .36E-08 

5.18E-09 

2.48E-09 

2.43E-10 

0. 

.20-. 25 

1 . 96E-08 

1. 18E-08 

5.83E-09 

2.26E-09 

9.51E-10 

.25-. 30 

2 .98E-08 

1. 57E-08 

7.70E-09 

3.07E-09 

1.18E-09 

. 30-.35 

3 . lOE-08 

2.29E-08 

1. 18E-08 

4.96E-09 

2.12E-09 

.35-. 40 

4.27E-08 

3.09E-08 

1.68E-08 

6.09E-09 

9.R7E-] 0 

.40-. 50 

2 .84E-08 

3. 08E-08 

2. 13F-08 

6.52E-09 

1 .81E-09 

.50-. 60 

8.52E-08 

8. 09E-08 

5.68E-08 

1 .93E-08 

9 .90F-09 

.60-. 70 

9.69E-08 

1.03E-07 

7.44E-08 

3.22E-08 

1.31E-08 

.70-. 80 

1 .26E-07 

1.36E-07 

1 . 16E-07 

4.97E-08 

1 .85E-08 

.80-. 90 

1 .42E-07 

1. 49E-07 

1. 29E-07 

6.35E-08 

3 .06E-08 

.90-1.0 

1 .65E-07 

1.66E-07 

1.46E-07 

7.95E-08 

3.82E-08 

1.0-1. 5 

2.05E-07 

2.30E-07 

2.00E-07 

1 .05E-07 

5.64E-08 

1 .5-2 .0 

2.79E-07 

2.96E-07 

2.72E-07 

1 .53E-07 

1 .05E-07 

2. 0-2. 5 

3. 16E-07 

3.60E-07 

3.26E-07 

2.21E-07 

1.13E-07 

2. 5-3.0 

3 .54E-07 

4.28E-07 

3.85E-07 

2.58E-07 

1 .60E-07 

3. 0-3. 5 

4.33E-07 

5. 15E-07 

4.74E-07 

3.12E-07 

1.99E-07 

3. 5-4.0 

4.25E-07 

5.08E-07 

5.01E-07 

3.56E-07 

2.09E-07 

4.0-4. 5 

4.88E-07 

6.07E-07 

5.71E-07 

4.09E-07 

2.65E-07 

4. 5-5.0 

5 .21E-07 

6.63E-07 

6.50E-07 

4.50E-07 

2.65E-0 7 

5. 0-6.0 

6.03E-07 

7.31E-07 

7.06E-07 

5.16E-07 

3.32E-07 

6. 0-7.0 

6.61E-07 

8.61E-07 

8.39F-07 

5.82E-07 

3.85E-07 

7. 0-8.0 

7.69E-07 

9.99E-07 

9.25E-07 

6.76E-07 

4.71E-07 

8. 0-9.0 

9.06E-07 

1.08E-06 

1.07E-06 

7.86E-07 

5.08E-07 

9.0-10.0 

9.42E-07 

1. 18E-06 

1.17E-06 

8.89E-07 

5.63E-07 








TABLE VI.- TOTAL EXTREMITY DOSE RESPONSE FUNCTIONS 
DUE TO INCIDENT PROTONS 


Energy 

20 g/cm2 

58 g/cm2 

100 g/cm2 

200 g/cm2 

300 g/cm2 

range, 






GeV 

Dose response function, — „ 

incident proton/cm^ 

. If^-. 15 

1 .29E-09 

7. 71E-10 

2.81E-10 

0. 

0. 

. 15-.20 

9.68E-08 

1.03E-09 

4.95E-10 

4.86E-11 

0. 

.20-.25 

1.58E-07 

4. 13E-09 

1. 16E-09 

4.50E-10 

1 .90E-10 

.25-. 30 

1 . 52F-07 

6.69F-09 

1. 54E-09 

6.12E-10 

2.35E-10 

.30- .35 

1 .44F-07 

5. lOE-08 

2.36E-09 

9.89E-10 

4.23E-10 

.35-. 40 

1 .29E-07 

8.49F-08 

6. 91E-09 

1.21E-09 

1 .97E-10 

.40-. 50 

7.41E-08 

4.90F-08 

1.59E-08 

1 .30E-09 

3 .62E-1 0 

.50-. 60 

1 .63F-07 

1.C7F-07 

6.30F-08 

3.84E-09 

1 .97E-09 

.60-. 70 

1.75E-07 

1. 27E-C7 

7.18E-08 

1 .43E-08 

2 .60E-09 

.70-. 80 

1.91E-07 

1. 35E-07 

8.97E-08 

2.13E-08 

5.88E-09 

.8 0- .<70 

1 .94E-C7 

1.41F-07 

9. 54E-08 

3.65E-08 

9.12E-09 

.90-1.0 

2 .12E-07 

1 .46E-07 

1.22E-07 

5.50E-08 

1 .47E-08 

1.0-1. 5 

2.23F-07 

1 .85E-07 

1.40E-07 

6.57E-08 

2.48E-08 

1.5- 2.0 

2.81F-07 

2.46E-07 

1. 75E-07 

9.34E-08 

4. 71E-08 

2.0-2 .5 

3.41E-07 

2.64E-07 

2.12E-07 

1.19E-07 

6.95E-08 

2. 5-3.0 

3 .48E-07 

3. 20E-07 

2.47E-07 

1 .53E-07 

7.61E-08 

3. 0-3. 5 

3 .94E-07 

3.73E-07 

3.30E-07 

1 .80E-07 

1 .llF-07 

3 .5-4.0 

4.D6E-07 

3.90E-07 

3. 58E-07 

2 .39E-07 

1 .43E-07 

4. 0-4. 5 

4 .87E-07 

4.57E-07 

3.83E-07 

2.61E-07 

1.20E-07 

4.5-5. 0 

5.09E-Q7 

5. 19F-07 

4.42E-07 

2.65E-07 

1 .60E-07 

5 .0-6.0 

5.27E-07 

5.41E-07 

4.62E-07 

3.40E-07 

1 .81E-07 

6. 0-7.0 

6 .49E-07 

6.37E-07 

5.57E-07 

3.75E-07 

2.13E-07 

7. 0-8.0 

6.73E-07 

7. 21F-07 

6. 80E-07 

4.31E-07 

2.86E-07 

8 .0-9.0 

7.31F-07 

8.21E-07 

6.76E-07 

4.83E-07 

2.94E-07 

9. 0-10. 0 

8.38E-07 1 8.79E-07 

8. 20E-07 

5.05E-07 

3 .39E-07 



rem 





incident proton/cm 2 

. 10-. 15 

6.45E-09 

3.87E-09 

1.41E-09 

0. 

0. 

.15-. 20 

2 .43E-07 

5. 18E-09 

2.48E-09 

2.43E-10 

0. 

. 20-. 25 

3 .41E-07 

1.61E-08 

5.83E-09 

2.26E-09 

9.51E-10 

.25-. 30 

2 .88E-07 

2.44E-08 

7.70E-09 

3.07E-09 

1.18E-09 

.30-. 35 

2 .71E-07 

1. 33F-07 

1.18E-08 

4.96E-09 

2.12E-09 

.35-. 40 

2 .45E-07 

1.88F-07 

2.55E-08 

6.09E-09 

9.87E-1 0 

.40-. 50 

1 .44E-07 

1. 13E-07 

4.39E-08 

6.52E-09 

1 .81E-09 

.50-. 60 

3 .45E-07 

2.49E-07 

1.59E-07 

1 .93E-08 

9.90E-09 

.60-. 70 

3. 70F-07 

3. OOE-07 

1.75E-07 

5.04E-08 

1 .31E-08 

,70-. 80 

4. 22E-07 

3.40E-07 

2.35E-07 

7.01E-08 

2.33E-08 

.80- -90 

4. 28E-07 

3.43E-07 

2.60E-07 

1 .lOE-07 

3.64E-08 

.90-1 .0 

4.89E-07 

3.69E-07 

3.20E-07 

1 .56E-07 

5.34E-08 

1 .0-1.5 

5 .27F-07 

4. 83E-07 

3.88E-07 

1.90E-07 

8.22E-08 

1 .5-2.0 

6.92E-07 

6.41E-07 

4. 86E-07 

2.68E-07 

1 .53E-07 

2. 0-2. 5 

8.50E-07 

7. 16E-07 

5.98E-07 

3.66E-07 

2.08E-07 

2 .5-3.0 

8.84E-C7 

8. 74E-07 

7.04E-07 

4.51E-07 

2.42E-07 

3. 0-3. 5 

1 .D2E-06 

1.02E-06 

9. 37E-07 

5.45E-07 

3.39E-07 

3. 5-4.0 

1 . 03E-06 

1.07E-06 

l.OOE-06 

6.89E-07 

4.18E-07 

4. 0-4. 5 

1.25E-06 

1. 25E-06 

1.08E-06 

7.58E-07 

3.87E-07 

4 .5-5.0 

1.29E-06 

1.41E-06 

1.25E-06 

7.92E-07 

4.77E-07 

5. 0-6.0 

1 . 37E-06 

1.49E-06 

1.32E-06 

9.88E-07 

5.61E-07 

6. 0-7.0 

1 .67E-06 

1. 76F-06 

1.60E-06 

1.08E-06 

6.44E-07 

7. 0-8.0 

1 .78E-06 

2.02E-06 

1.9CE-06 

1 .25E-06 

8.53E-07 

8. 0-9.0 

1 ,97E-06 

2. 28E-06 

1.98E-06 

1.41E-06 

8.84E-07 

9.0-10.0 

2 . 19E-06 

2.43E-06 

2. 31E-06 

1.52E-06 

1 .OlE-06 
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Figure 1.- Proton cutoff energy for vertical and western horizon 
as a function of geomagnetic latitude. 



Transmission fraction 






injection angle, a, deg 

(b-3) Accelerating voltage, 800 volts; injection angle, -o. (b-4) Accelerating voltage, IZOOvoIts; injection angle, -a. 

(b) Atidplane injection radius, 4 centimeters. Curve sho^vn is transmission fraction calculated from beam divergence 
measurements for loss cone angles of 42® and 43®. 
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Figure 16. - Concluded. 
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experimental data and the predicted results. 

In figure 16(b), transmission results are presented for electrons injected at the 
midplane at a radial distance of 4. 0 centimeters from the axis and for a range of condi- 
tions. The mirror ratio (taking into accomit the curvature of the magnetic field 
lines) at a radial distance of 4,0 centimeters is 0.446. The loss cone angle associated 
with this mirror ratio is 42° (eqs. (6)). 

In figure 16(b-l) and (b-2) the transmission results are shown for positive injection 
angles. The curve calculated from the beam divergence and a loss cone angle of 42° is 
also plotted on the figures for comparison. Because electrons injected off-axis at posi- 
tive injection angles gyrate above the injection point, the mirror ratio for this region 
slightly less than 0. 446 resulting in a loss cone angle slightly less than 42°. The ex- 
perimental results again agree closely with those predicted from beam divergence mea- 
surements. 

In figure 16(b-3) and (b-4) the data are presented for negative injection angles, 
where the electrons gyrate below the injection point. In this region the mirror ratio ap- 
proaches the axial mirror ratio of 0. 465 which corresponds to a loss cone angle of 43°. 
Calculated results are plotted for loss cone angles of 42° and 43°. From theoretical 
considerations, the data points would be expected to fall between these transmission 
curves. The actual data points, however, are distributed about the transmission curve 
for 43°. This slight shift may be data scatter or may be the result of nonadiabatic be- 
havior. For the range of e (0.029 to 0. 162) covered by the data of figure 16(b), the loss 
cone angle did not increase (within the 2° spread of the data). 

The transmission results presented in figure 16(c) were obtained with an electron 
gun positioned on axis (R = 0) but slightly behind the midplane (z = -0. 5 cm). These re- 
sults were used to identify any change in the transmission characteristics between a 
particle crossing the midplane and a particle injected at the midplane. A change in the 
transmission characteristics might be expected because the largest change in the mag- 
nitude of the magnetic moment occurs at the midplane. The electron gun used for these 
measurements was a triode gun with an oxide coated cathode. An anode aperture of 
0. 15 centimeter resulted in a beam divergence angle as large as ±15°. 

For the tests plotted in figure 16(c-l) and (c-2) the magnetic field was held constant 
while the beam accelerating voltage was varied. The beam angular distribution in- 
creased with accelerating voltage as evidenced by the changing slope of the transmission 
measurements. The slope of the data is constant for each accelerating voltage though, 
and passes through the 50 percent transmission point at a loss cone angle of 43°. The 
range of e for figure 16(c-l) extended from 0. 04 to 0. 07, and for figure 16(c-2) from 
0.031 to 0.055. Within this limited range of €, the results indicate that the loss cone 
angle did not increase as a result of injection slightly behind the midplane. 

In summary, the experimental range of e investigated in the course of this effort 
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was from 0.027 to 0.22. Results of the transmission measurements (fig. 16) indicate 
that the loss cone did not increase (within the 2° spread of the experimental data). For 
single interaction with the magnetic mirror region, particle behavior was adequately 
predicted by adiabatic theory. Values of e considerably larger than 0.22 were not in- 
vestigated with the experimental apparatus described herein because the required low 
magnetic field strengths (resulting in large electron gyro radius) would lead to electrons 
striking the magnet coil. This was found by test to occur at an e = 0. 35. Values of e 
below 0.02 required magnetic field strengths that distorted the distribution of the elec- 
tron beam. 

Previous results (ref. 13) for single interaction with a magnetic mirror using much 
larger magnetic fields indicated that a transition to nonadiabatic behavior would be ex- 
pected for e = 0. 046 at R = 0 and for e = 0.041 at R = 4.0 centimeters for a mirror 
ratio of 0. 465. In the range of values tested and reported herein, which extends both 
above and below the indicated adiabatic transition value of e, no evidence of an increase 
in the loss cone was found. 


CONCLUDING REAAARKS 

Nonrelativistic single particle behavior was experimentally investigated in an axi- 
symmetric mirror magnetic field. Measurements were made on an electron beam for 
a single interaction with the mirror region. Within the 2*^ spread of data of this investi- 
gation, the loss cone angle for a single interaction of electrons with a mirror field was 
not observed to increase (i. e. , to become noticably nonadiabatic) for a range of the adia- 
batic parameter e from 0.027 to 0.22. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 25, 1970, 

120-26. 
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’ incident proton/cm 



E,GeV 

Figure 5-- Proton extremity dose response functions due to incident protons 
at various altitudes. (Note the Bragg peak at 20 g/cm2.) 




incident proton/cm 



E.GeV 

Figure 6.- Neutron extremity dose response functions due to incident protons 

at various altitudes. 
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Altitude, g/ cm 

Figure 11.- Upper and lower limits of dose rate in 
extremities for the prompt spectrum of the 
Fehruary 195^ solar event . 
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Figure 12 .- Dose rate averaged over the whole "body 
for the solar event of irovemher 12, i960. 
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Figure 13 •- Dose rate in extremities for the 
solar event of November 12, i960. 
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